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Abstract 

The behavior of a massive star during its final cata- 
strophic stages of evolution has been investigated theoreti- 
cally, with particular emphasis upon the effect of electron- 
type neutrino interactions. The methods of numerical hydro- 
dynamics, with coupled energy transfer in the diffusion 
approximation, were used. Gravitational collapse initiated 
by electron-capture and by thermal disintegration of nuclei 
in the stellar center is examined, and the subsequent be- 
havior does not depend sensitively upon which process causes 
the collapse. 

As the density and temperature of the collapsing stellar 
core increase, the material becomes opaque to electron- type 
neutrinos and energy is transferred by these neutrinos to 
regions of the star less tightly bound by gravity. Ejection 
of the outer layers of the star can result. This phenomena 
has been identified with supernovae. 

Uncertainty concerning the equation of state of a hot, 
dense nucleon gas causes uncertainty in the temperature of 
the collapsing matter. This affects the rate of energy trans- 
fer by electron- type neutrinos and the rate of energy lost 
to the star by muon-type neutrinos. 


f 


II 

The effects of general relativity do not appear to become 
important in the core until after the ejection of the outer 
layers . 
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Introduction 

The behavior of a massive star during its final catastrophic 
evolution has been investigated theoretically, with particular 
emphasis upon the effect of electron-type neutrino interactions. 
Colgate and White (1964) have suggested that the gravitational 
collapse of such a star may be partially reversed by a combination 
of shock phenomena and energy transfer by neutrino diffusion from 
a hot, ultra-dense core. The resulting ejection of hot matter 
has been identified with supernovae. The Von Neuman-Richtmeyer 
pseudo-viscosity method of numerical hydrodynamics, coupled with 
energy transfer in the diffusion approximation, has been used to 
investigate this hypothesis. 

In section I the physical processes involved in the collapse, 
and the gravitational stability of a massive star are discussed. 

In particular, the problems of constructing an equation of state 
and determining the energy transfer by neutrinos under the extreme 
densities and temperatures to be encountered are considered. 
Following this, section II develops the initial models and presents 
their subsequent histories. Three methods of treating neutrino 
energy transfer - (1) no energy loss, (2) energy loss by electron 
pair-annihilation and plasmon decay neutrinos, and (3) thermal 
diffusion of neutrinos - are presented, and the results contrasted. 
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The effect of initial structure upon subsequent history is 
examined. The behavior of these models is compared with the work 
of Colgate and White (1964). 

Section III contains an examination of the implications of 
the calculations reported in II, and section IV is a critique of 
the methods employed in the calculations. Details of the numerical 
techniques of hydrodynamics and energy transfer are presented in 
an appendix. 
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1. Physical Processes 

In order to maintain its luminosity a star may derive energy 
from two sources: thermonuclear reactions and gravitational 

contraction. In the latter case the star converts its gravitational 
potential energy into kinetic energy of gas particles on such a 
slow time scale that hydrostatic equilibrium is approximately 
satisfied. No nuclear fuel can last indefinitely, so one expects 
that eventually the star will contract to higher densities. At 
these higher densities the Pauli exclusion principle can become 
operative and contribute to the pressure. However Chandrasekhar 
(1939) has shown that for electrons the maximum mass of a body 
supporting itself by degeneracy pressure is less than 1.5 solar 
masses.* If the mass of a star is less than the Chandrasekhar 
limit, it may radiate away its remaining thermal energy and settle 
down as a white dwarf. For more massive stars the situation is 
not so simple. 

A massive star spends most of its life burning hydrogen and 

helium.** Helium-burning produces oxygen and perhaps some carbon 

(Deinzer (1964) ). Fowler and Hoyle (1964) have discussed in 

detail the nuclear reactions occurring in subsequent evolution. 

* For pure hydrogen the limit is higher, (see Chandrasekhar (1939), 
p. 423) but a pure hydrogen star is unrealistic when the Fermi 
energy of electrons is higher than the beta-decay energy of 
the neutron. 

** See Hayashi (1962), Hofmeister (1964), and Stothers (1965). 
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After carbon-burning (temperature~8xl0® °K) neutrino production 
by plasmon decay and electron-positron pair-annihilation robs the 
star of significant amounts of energy. This speeds the evolution 
of the star. Chiu (1964) has calculated some models of pre- 
supernova stars including this effect. 

There are at least two mechanisms by which the star can be 
robbed of internal energy faster than it can replace the lost energy 
by quasi-static gravitational contraction. They are thermal 
disintegration of nuclei and electron-capture. When these 
mechanisms operate, the star will collapse, falling almost freely 
in its own gravitational field. Which process will actually trigger 
the collapse depends on the details of pre-implosion evolution. 

A. Thermal Disintegration of Nuclei 

For temperatures greater than T = 4x10^ °K and densities the 
order of or greater than p = 10^ gm/cm^ a wide variety of nuclear 
reactions can occur. A calculation of these rates requires an 
accurate knowledge of the initial nuclear composition of the matter, 
a large collection of nuclear parameters,* and considerable effort. 
It is beyond the scope of this work to justify a pre-supernova 
model involving such a calculation. 


* See Truran, Hansen, Cameron, and Gilbert (1965), for instance. 
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If the reactions are fast enough, the problem may be treated 
by the methods of statistical equilibrium. Insofar as the 
equation of state is concerned, the only changes of interest are 
those which are strongly endoergic or exoergic. The photo- 
disintegration of ^ 8 Fe matter in the implosion has been discussed 
by Fowler and Hoyle (1964) . They find that for temperatures 

T 7. 10 9 °K that the photodisintegration time is t ~ 10 -6 second 

which is no larger than the most restrictive hydrodynanu^ time 

scale. Iloyle and Fowler note that the equilibrium composition for 

56 4 

Fe matter changes to essentially pure He in a region of width 

A T ~ lxlO 9 °K 

about the density-temperature curve corresponding to an equilibrium 

. 56 4 

concentration of half Fe, and half He and neutrons, i.e., the curve 

39.17 

log p = 11.62 + 1.5 log T^ - Tg 
where T^ means temperature in units of one billion degrees, 
logarithms are to the base 10 and density p is in gm/cm 3 . This 
may be approximated by the expression 

T 9 /6.0 = (p/1.82 x 10 6 gm/cm 3 ) 0 * 081 . 

For somewhat higher temperatures photodisintegration of alpha 
particles is expected. 

4 He - 2p + 2n 


4 
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In this case the transition region is approximated by 
T/(12xl0 9 °K) = (p/(10 8 gm/cm 3 ) )°- 13 
Now the energy required to produce the reaction 

56 Fe - 13 4 He + 4n 
is 

Q(Fe,a) = -2.1x10^-® erg/gm 

and for 

4 He - 2p + 2n 
is 

Q(a»np) = -6.8x 10 13 erg/gm 

These values will be roughly correct even if the composition is 
not ^Fe, but some other stable nuclide Subsequent results will 
not depend sensitively upon this choice. 

B. Klee tr on-capture 

Consider a nucleur Z,A) which is stable aga’ist beti.-decay 
or earth. It may be a product of some react ic" , epresented by 

(1) (Z-l, A) r*. is, A) + e” + 

Under stellar conditions of extreme density an endothermic 
reaction of the form 

(2) (Z, A) + e“ - (Z-l, A) + v 

can occur in which the terrestrially stable nucleus (Z,A) is 
induced to captui ) a continuum electron from the surrounding plasma. 
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Bahcall (1964c) has investigated this process and finds that 
for allowed decays, the half-life for the process (2) in a 
stellar interior is related to the half-life for the decay 
(1) on earth by 


^ T 1/2 ^ star (fT l/2 ) lab K 


-1 


where f is the usual function used 


-he comparative half 


life, i.e.. 


f (±Z,W) 


rmax , 2 . „ TT , 

! dp p F(+Z,W) 
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and p is electron momentum, W is relativistic electron energy, 
q is the neutrino momentum and F(+Z,W) is the Fermi function. 


Now 

, 2 2 / w . (W-|i)/kT) 

K ** j dp p q F(Z,W)/(l+e 

p 

threshold 

(where n is the relativistic chemical potential) i*; the appro- 
priate generalization of f(+Z,W) for reaction (2). The range 
of integration now extends from the electron momentum corre- 
sponding to the threshold energy of (2) to all higher enernies. 
Tho correct weighing factor for a Fermi-Dirac distribution of 
electrons , 

(1 + exp [ (W-p)AT] '“ 1 , 

is included. The integral K will be large compared to f(+Z,W) 

when the electron distribution is such that energy levels with 

W > W .. . . _ 

threshold 

for the capture reaction, are well populated. This can occur in 
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the nondegenerate case when 

(3) kT > W 

threshold 

and in the degenerate case when 

(4) W > W 

Fermi threshold 

Taking ^ to t ^ ie or <3er of nucleon binding energy in 

the nuclear potential should give an estimate of the thermodynamic 
conditions under which induced electron capture will begin to 
occur. If we take the threshold for the electron capture to 
be 


W ~ 8 MeV, 

threshold 

then for a nondegenerate gas, the condition (3) implies that the 
temperature is 

T > 10 11 °K, 


but thermal disintegration discussed in the previous section will 
have already disrupted the nuclei at much lower temperatures. Using 

(4) gives a condition* on the density p 

p 9 3 

“ £ 10 gm/cm . 

M e 


* The quantity |j has its usual astrophysical definition of average 
atomic weight per free particle ( ye _ + for a completely 

ionized gas) . 
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Beyond this density nuclei will capture high energy electrons 

4 

from the top of the Fermi sea. The alpha particle. He, will 

be extrenu.ly ..esistent to electron capture because there is 

no stable nucleus with Z = 1, A = 4, whereas the a. particle 

4 

is tightly bound. The threshold for electron capture on He 
will be of the order of the energy needed to disintegrate 
the nucleons, i.e., about 30 MeV for complete disruption, 

but the reaction of most importance is probably 

4 - 3 

He + e - H + ntv 

e 

which requires about 21 MeV. A value for the electron Fermi 
energy of 30 MeV corresponds to a density of the order of 


P 0 lrt ll . 3 

— ~ 2 x 10 gm/cm 

^e 

As the density rises, so that the Fermi energy becomes greater 
than the threshold energy, the continuum electron capture rates 
increase until an assembly composed predominantly of neutrons is 
formed . 


C. Stability Against Continued Implosion 
Once the implosion begins it cannot be stopped until the 
pressure is again large enough to provide mechanical support for 
the configuration. Chandrasekhar (1939) has shown that the 
boundary for mechanical stability of a self-gravitating mass is 
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the Y = 4/3 adiabat. That is, if upon compression the change 
in pressure and density of the material can be represented by 

Y = d(log P)/d(log p ) 

then for y > 4/3 the material is stable, but for Y < 4/3 the 
material is unstable toward continued contraction.* As particles 
become relativistic the relation between energy and momentum 
changes from 

e = sl 
2m 


to 


e = pc 

in the extreme relativistic case. The corresponding relations 
for a gas of such particles are 

nonrel . rel . 


E = 3 PV 
2 


E = 3 PV 


Y = 5/3 


Y = 4/3 


♦Actually the criteria are somewhat more complicated, dealing with 
pressure averages of Y.. F. Dyson, "Hydrostatic Instability of 
a Star," unpublished. 
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where E is the energy density per unit mass, P the pressure and 
V = 1/p the specific volume. Thus, as the temperature (or the 
Fermi energy) increases we expect the effective adiabatic exponent 
of the gas,Y , to approach four-thirds. At what temperature this 
occurs depends upon the rest mass of the gas particles (photons, 
having zero rest mass are always relativistic, electrons become 
relativistic for T ^ 6x10^ °K, while nucleons require T ~ 1.2 x 10-^ 
°K) . Large amounts of energy are removed from the star by neutrino 
emission. Both the processes of electron-capture and thermal 
disintegration of nuclei require large amounts of energy to 
proceed; in both cases this energy is at least of the order of 
nuclear binding energy. In view of these large energy requirements, 
the material is expected to become degenerate even if it was not 
originally so. 

The problem of stability against continued collapse is then 
reduced to whether a cold, dense neutron gas can give a pressure 
contribution which increases with density faster than Y= 4/3. 

From investigation of the properties of nuclear matter such a 
contribution is found. 

D. The Equation of State 

At densities of the order of or less than nuclear densities 
(p < 3 x 10^^ grn/cm^) the attractive nuclear potential lowers 
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the pressure below that expected for a degenerate, noninteracting 
gas of fermions. The nuclear potential becomes strongly repulsive 
at higher densities, and raises the pressure above that expected 
for a noninteracting gas, but the exact details of nuclear potentials 
in this range (greater than nuclear density) is not well known. 

From the several forms of the nuclear potential discussed by 
Tsuruta (1964), it appears that there will be a pressure term of 
the form 



wr*ere 

Y> 2 

These results for the equation of state are based on the assumption 
that the nucleons may represent a noninteracting, degenerate gas 
of Fermi particles in a common potential well. 

Bahcall and Wolf (1965) have attempted to determine the effect 
of nucleon-nucleon interactions more accurately by using the 
"independent-pair" model of Gomes, Walecka and Weisskopf (1958). 
This technique is valid only if the nucleons are highly degenerate. 
Unfortunately it is necessary to know the equation of state for 
nondegenerate and semi-degenerate nucleon matter. In view of the 
uncertainties involved in any nuclear equation of state and the 
numerical limitations of this investigation, and extremely simple 
form for nuclear pressure at high density and temperature was 
chosen: The nucleons were assumed to be a gas of noninteracting. 
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free Fermi particles. This reproduces the correct general 
character in the limits of complete degeneracy and of high 
temperature , low density. In order to avoid excessive use 
of computer time, the equation of state was constructed 
from a composite of analytic terms. 

In addition to the nucleon pressure terms discussed 
above, black-body radiation pressure and electron pressure 
(including relativistic degeneracy) were taken into account. 
Thus the approximate expression for the pressure is 

p - R dr) T + K nQ~) 

n n 


+R 





where R is the gas constant, the constants K and K are 

e n 


** 1.201 x 10 ^ dynes cm ^ 

= 5.226 x 10^ dynes cm ^ 

-3 

if the density p has units gm cm , T is the temperature, 

and the number density of the i — type of particle (n * neutron, 

e * electron) is 


where N is Avagadro's number. The number density of neutrons 

cl 

is negligible before electron capture occurs. Electron-pair 
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creation does not affect the equation of state for large electron 
number density. The energy density corresponding to this pressure 
has the simple form 




(v-D 


where y = 4/3 for the relativistic particles and Y - 5/3 for 
nonrelativistic ones. 

Electron capture reactions were accounted for as follows: 
Since the expression 


£_ - i 

He V U e 


n 


is proportional to the number density of electrons, smoothly 
changing u provides a convenient way to reduce the electron 
pressure of the system. When the relativistic electron Fermi 
energy E^, 






2 

m c 
e 

where p c is density in units of 10® gm/cm^, reaches a given 

O 

level, the number density of electrons is held constant until the 
degenerate nucleon pressure becomes more important. That this 
agrees with other estimates may be seen in Figure (1) . Two 
parameters are involved: the Fermi energy at which captures are 
supposed to occur, and the factor by which the electron number 
density is decreased. 


:mperature equation of state 



Density (g/cm 
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Actually the onset of electron capture reactions is a gradual 

9 3 

process, beginning at densities as low as 10 gm/cm . Figure (1) 

shows that the Harrison-Wheeler* equation of state gives a 

slightly lower pressure than expected from a relativistic electron 

gas at densities from about 10^ to lO^ gm/cm^ because of electron 

capture. Tsuruta (1964) finds that the equilibrium abundance 

5 6 

peak of zero temperature matter shifts from Fe at low density 

7 3 . ... 

(''■’10 gm/cm ) to very neutron-rich heavy nuclei at higher density 

11 

(>10 gm/cm ) due to electron capture. The effect of this is 
to initiate a contraction in a pre-supernova star at much lower 
densities than might be expected otherwise. The assumption of a 
sharp electron capture threshold is unrealistic, but is probably 
a minor source of error. In the calculations to be reported, 
electron capture was assumed to occur rapidly when the Fermi 
energy was slightly above the 21 MeV or so necessary to drive 
the reaction 

4 - 3 

He + e J H + n 

that is, for p ~ 10^ gnv/cm^ . 


* Harrison, B. K. , et al, 1965, chapter 10. 
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As Figure (1) shows, a zero-temperature body undergoing 

quasi-static gravitational contraction will encounter a pressure 

11 3 

deficit at densities above p ~ 10 gm/cm . The equation of ' 

state does not compensate for this pressure loss relative to a 
Y = 4/3 adiabat until the density is about 10^ gm/cm^. This is 
accomplished here by approximating the complicated equation of 
state for nuclear matter by that of a noninteracting non-relativistic 
Fermi gas. In the region where the pressure deficit occurs, 

10^ < p < 10^ gm/cm^, the detailed nature of the equation 

of state is relatively unimportant because the supernova core is 
falling in freely to higher densities and is not affected by the 
nucleon pressv *e contribution. In this region Tsuruta’s equation 
of state gives pressures between those of the Harrison-Wheeler 
and the Salpeter equation of state, except when a nuclear hard-core 
term begins to dominate the equation of state. It should be 
emphasized that the zero-temperature equation of state is only 
a convenient limiting case, and that finite temperature effects 
are important in super novae collapse. 

E. Neutrinos and Energy Transfer 
The effectiveness of a given mechanism for energy transfer 
depends on the rate at which energy can be put into the given mode, 
and on the speed with which the given mode moves this energy. 

Energy transfer in stars is generally accomplished by photon 
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diffusion, or in some cases by convective mass motion or conduction 

/ 

by degenerate electrons (see Schwarzschild , 1958). Using the 
Thomson cross-section for the electron, the photon mean free 


path is roughly 


1 photon - 1_ 


1 cm. 


N ' . v, 

e th 


where N e is the electron number density per unit volume, a th is the 

3 

Thomson cross-section, and p is the density in gm/cm . For the 

/* O 

core of a star with a central density of, say, 1C D gm/cm , the time 

for a photon to diffuse through even ten kilometers of matter is 

2 

T ph.dif. ~ AR ~ 10 years. 


The universal theory of weak interactions* of Feynman and Gell-Mann 
(1958) predicts a large number of processes that result in the 
formation of neutrino-antineutrino pairs. The emission of a 
neutrino pair is much less probable than the emission of a photon 
so that the process is not generally observable in the laboratory. 
Once formed, however, the neutrino pair is virtually certain to 
escape from a normal star ( p centra j_ « 10 gm/cm ) . For temper- 
atures less than several billion degrees, the cross-section for 
neutrinos and antineutrinos is roughly 

~ 10" 44 cm 2 . 


so that the mean free path is 


neutrino ~ 10 20 
P 


*Or ary theory which predicts coupling of terms of the form (ve) (ev). 


18 


which is roughly 100 light years for the density of the sun. 

Thus, excluding extreme conditions, the production of neutrinos 
acts as an instantaneous local energy sink for the star. 

1. Energy Loss by Neutrino Escape. 

Chiu (1961) has calculated the energy loss rate due to the 
process 

e + + e” v + v 

for stellar material in the temperature range (0.5 to 10) x 10^ °K 
and densities (0 to 10^ gm/cm^) . In much of this range the electrons 
are partially degenerate and numerical eveluation of integrals was 
necessary. Analytic forms for limiting cases hava been presented 
by Chiu and Stabler* (1961) . 

For 

mc^ << kT 

and 


E . £ kT 

Fermi 


the energy loss rate is 


Q = 4.3 x 10 


15 


(* 9 ) 


ergs/gm/sec 


where p is in gm/cm J and T 9 is the temperature in units of 10 & °K. 


* See Ritus (1962) for a numerical correction of the photoneutrino 


rates in this paper 
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This approximation is inaccurate when temperatures are about 
9 o 

3 X 10 K or below, but the energy loss rates are then too small 
to affect the hydrodynamic calculation. The approximation is also 
invalid when degeneracy is pronounced, but then the plasmon-neutr ino 
loss rate is larger, so that the numerical error is negligible. 

The production of neutrino pairs by coherent electron ex- 
citations (transverse plasmons) in a hot, partially degenerate 
relativistic plasma has been calculated* by Adams, Rudermar. and 
Woo (1963) , and extended by Inman and Ruderman (1964) . Neutrino- 
pair emission by collective electron modes, especially transverse 
plasma excitations, is found to be the main mechanism for neutrino 
radiation by a dense plasma when electron-positron production is 
small either because the temperature is too low or degeneracy 
suppresses it. Chiu (1964) gives an analytical approximation 
for the plasma neutrino process: 

Q * - 1.1 (T 9 ) 3 p erg/gm/sec. 

where the usual notation is employed. This is valid for 

x <: 1 


where x is given by 

2/3 ) -1/4 1/2. 

x = 0.237(1 + 0.6413( P6 ) (p & ) 

if 


P 6 » 1 

* Zaidi has recently indicated that this rate is too large by a 
factor of 4 (to be published) . 
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r *0 

The density p ^ has units of 10° gm/cm . Actually numerical 

evaluation of the pair annihilation, plasma and photoneutrino 

rates by Hansen (1964) indicates that this approximation is 

11 1 

reasonably good (factor of 2) for densities as high as 10 gm/cm 

or so if the temperature is T ^ 10 x °K. Again, errors in small 

energy-loss rates are negligible from a hydrodynamic point of view. 

For temperatures much higher than this the pair annihilation rate 

is dominant, so that the approximation is reasonable in the region 

in which it is the primary energy loss mechanism. 

The plasmon-decay and the pair-annihilation energy loss rates 

are of interest for two reasons. First, they may be important 

in cooling shock-ejected matter whose density is less than 
11 3 

p ~ 10 gm/cm . Also, for higher densities other neutrino 
processes will probably dominate, but these processes nevertheless 
provide a convenient estimate of the minimum possible energy loss 
rate which is not plagued by the uncertainties* in the effect of 
strong interactions. 

2. Neutrino Opacity. 

The interactions of neutrinos ard antineutrinos in dense matter 
have been discussed by Bahcall (1964a), Bahcall and Frautschi (1964b) 
and Euwema (1964) . Euwema calculates the inhibiting factor for 
effects of the exclusion principle on neutrino absorption. 

* See Bahcall and Wolf (1965) for example. 
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Negligible temperature and completely noninteracting particles 
which filled all states below the Fermi level and none above, were 
assumed. Bahcall and Frautschi have considered neutrino-lepton 
scattering and neutrino-nucleon interactions generally. In 
particular, Bahcal?. (1964a) has suggested that neutrino-electron 
scattering 

/ / 

v 9 + e” - ^8 + e" 

/ / 

+ e “* Vg + e” 

is the most important neutrino process for energy deposition in the 
supernova model of Colgate and White (1964) . 

a) Nondegenerate gas . 

For a nondegenerate gas of electrons, the total cross-section 
averaged over the initial electron distribution is 

3 -K 3 -1 3 r (W-nJATl.i 

<<Vp - [ (“> 1 1 J a 3 p Li + e - 

where a(p a ,uu a ) for neutrino-electron scattering is 

2 

° (Prt' ^a) = o o lPn' w a) /(! + 2 Pa ,l %) 
and for antineutrino electron scattering is 

a (Pa'%) = cr [1 - (1 + 2 Pa* UJ a)” 3 ] 

° v 6 

where n is the electron number density per unit volume, p is the 
e 

electron momentum, W the total electron energy, p a the dimensionless 
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four-momentum of the electron, ^ the corresponding four-momentum 

for the neutrino, and (j the chemical potential of the electrons. 

The constant a is 
o 


a =4 


o - (— )" 4 ( -%) 

tt me' n\ cr' 
e e 


-44 2 

^ .1.7 x 10 cm^ # 

2 

In the case kT » n» e c , the thermal motion of the electrons produces 
a large center-of-mass energy and hence causes the cross-section 
to exceed greatly the cross-sections for electrons at rest. In 


this limit Bahcall finds, for neutrino energies E*. » m c , that 

v e 


( a )uj,T - 3.2 (kT/m c ) 

e 2 


for the neutrino, and for the antineutrino 

a tu ct iw 

O' _ o 


Bahcall gives the following form for the neutrino energy loss 


per scatter 


(W-U)') 


av. 


UL) 


* 2 [l - ] 


cum c 
e 

where u) is the dimensionless neutrino energy and the prime refers 

to the quantity after collision. This approximation requires 

ou » 1 and kT » m c 2 

e 

For sufficiently high electron temperatures, the neutrinos gain 
more energy per collision on the average than they lose. 
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b) Degenerate Gas . 

In a degenerate electron gas both the initial electron 
distribution and the prior occupation of final electron states 
must be considered. Bahcall (1964a) gives the general expres- 
sion for neutrino-electron scattering as 


a 

U) 


= CT or4n 3 n s ft V^-T 2 

4ttL 6 C 2 ) J 

m c 
e 



P 

Tv 


dV 


S(e)Cl - S(€ / )](P C [ a (4) (P a +w a -p ( J-«» ( ^ ) 1 

where e is the dimensionless electron total energy and the 
other notation is as before, except for 

S (x) s [l + exp ([x-|i]AT)j 

where p is the electron chemical potential. For a completely 
deg merate gas , 


S( e ) = 1 for e < E f 
= 0 for e > 

where E_ is the total Fermi energy. Bahcall estimates 

1 2 

O uu U)«E 

_ .... / _ . « f o f 

a ( w ) 1 

*■ a E ui uu»E _ 
of f 

and for antineutrino scattering the results are multiplied 
by 1/3. 

The results of Bahcall and Frautschi (1964b) indicate 


that for 


n > 10 " 2 n 

e n 


a 


m 9 



i 
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that is, for the electron number density more than 1 % of the 
nucleon number density, neutrino-nucleon scattering is less than 
neutrino-electron scattering, and will be neglected. 

The situation is not clear concerning neutrino absorption by 
nucleons. In the ''low” density domain (p < 3x10 gm/cm ) nucleons 
not bound in nuclei will display the cross-section 

O ((D) „ 0 Q U) 2 

2 

for neutrinos with energy much larger than me ~ 0.5 MeV. For 
lower energies the behavior is more complicated, with the reaction 

v + p - n + e + 
having a threshold while the reaction 

v + n p + e~ 

does not. For densities greater than or the order of nuclear 
density ( p ;> 1x10 gm/cm J ) the previously mentioned uncertainty 
with the equation of state may affect the neutrino absorption 
cross-section by changing the threshold energy or by reducing the 
phase space available to nucleons in the exit channel. In any case, 
the extinction cross-section for neutrinos is at least as large as 
that predicted by neutrino-electron scattering alone, and may be 
larger. If the number density of electrons becomes considerably 
less than that of nucleons, then the neutrino opacity of the material 
cannot be represented, even approximately, without knowing the 
neutrino-nucleon interaction cross-sections. Thus the equation of 
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state problem is again encountered. The opacities used in the 
calculations presented here are due to neutrino-electron scattering. 

3. Energy Transfer by Neutrinos. 

In this investigation neutrinos were assumed to transfer energy 
by thermal diffusion. Although tested numerical techniques for 
time-dependent transfer problems are available (Richtmyer, 1957) 
it appears that any reasonably accurate transport treatment requires 
too much machine time for an exploratory calculation. In addition, 
other approximations, such as neglect of general relativity, make 
a detailed study of the coupled problem inappropriate at this time. 
Consequently energy transfer by neutrinos was treated in the thermal 
diffusion approximation* with the hope that the accuracy would be 
on the level attained by other aspects of the calculation. 

The thermal diffusion approximation assumes that the diffusing 
energy carriers (usually photons, but in this case electron-type 
neutrinos and their anti-particles) are in thermal equilibrium with 
the medium through which they move. The anisotropy which drives 
the energy transfer is assumed to be a perturbation on a generally 
isotropic distribution of carriers. The assumption of thermal 
equilibrium avoids detailed consideration of the processes of 
neutrino formation, a neglect which greatly simplifies the problem. 

* The approach in this section follows that of Frank-Kamenetskii 


(1962) for photons. 
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The first law of thermodynamics can now be expressed as 
dE = (s - ) dt - PdV 

where E is energy density, s the rate at which thermal energy is 
added to the medium, P is the pressure, and V is the specific volume. 
Assuming spherical geometry, the mass M satisfies the equation 

dM = 4 n p r^dr 

where p = 1/V is the density and r is the radius of the spherical 

element under consideration. The luminosity L is given by the 

usual expression 

L * (4rrr^) ^ ac d (T 4 ) 

3K dM 

where a is the radiation constant, c the velocity of light, T the 
temperature and K the Rosseland mean opacity 

K = 1_ 

pZ 

where is the corresponding mean free path. Notice that 
for simplicity we have assumed that the neutrino and anti- 
neutrinos may be described by a single "black-body" Fermi 
gas, that is, the chemical potentials of the neutrinos and 
antineutrinos are zero. Using Bahcall's (1964) limiting 
forms for the electron-neutrino (antineutrino) cross-sections, 
the Rosseland mean opacities can be evaluated analytically* 

* L.D. Landau and E.M. Lifshitz (1958) Statistical Physics. 
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for a "black-body" Fermi distribution. However it is 
sufficiently accurate to replace the neutrino energy by 
an average value 

< e > ~ 3kT 
v 

so that 

pK = N a (3kT) 2 
w e o 

where the quantities have been defined previously. 

It is interesting to note that the integrated neutrino- 
electron scattering cross-section for a degenerate electron 
gas is nonzero even if the neutrino energy u> is less than 
the electron Fermi energy e^, 

2 

< a > ~ a Q iu for (!)<<€£ 

This means that there is no completely transparent window 

2 

even for low energy neutrinos (m c « u) « c ) . Very low 

2 

energy neutrinos (uu < m e c ) are expected to transfer little 
energy on short time scales. 

The approximation of thermally diffusing neutrinos 
will be valid only if the neutrino mean free path is shorter 
than the distance in which the macroscopic variables change. 
That this condition may be satisfied can be seen as follows. 
For a density greater than 


p ~ 10^ gm/cra 2 
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the electron Fermi energy never falls below about 30 MeV, so 

that the electron number density is at least N ~ 10^5 cnT^. 

e 

Detailed numerical calculations show that the macroscopic 

variables change little over distances of the order of 

6 

A X < 2x10 cm. 

If the mean free path for a neutrino must satisfy the relation 

* £ A X 

then the average neutrino energy must be 

e £ (A X N e Co)’* 5 ~ 8 MeV 
If a thermal distribution is assumed, e ~ 3kT, so that the 
temperature must be greater than 

T ~ 36 x 10 9 °K. 

Temperatures far in excess of this are encountered. Eventually 
this condition breaks down at lower densities, and a "luminous 
surface" for neutrinos is formed, beyond which the neutrinos 
almost certainly escape the star without interaction. In this 
region the energy deposited by the incident neutrino flux and 
that lost by neutrinos escaping were taken into account in 
determining the boundary condition. Because of the rapid 
transition from neutrino-opaque to neutrino- transparent condition, 
the calculation appears to be insensitive to exact form used. 
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III. Hydrodynamic Calculations of Stellar Collapse. 

A. Initial Models. 

Models for stars just prior to supernova implosion have 
been suggested by Fowler and Hoyle, (1964) , and by Chiu (1964) . 
It appears that the essential differences are: (1) Chiu's 
model is much more centrally condensed, with a higher central 
density, and (2) because the central temperatures are roughly 
the same, the center of Chiu's model lies on a lower adiabat. 

The model suggested by Chiu has electron degeneracy in the 0013 
while that of Fowler and Hoyle is nondegenerate, at least 
until endoergic nuclear reactions become significant. Figures 
(2) and (3) illustrate the two approximate models chosen to 
reproduce these characteristics. An n = 3 polytrope, that is, 
a gravitating gas sphere in hydrostatic equilibrium for which 
the pressure and density are related by 

P ~ p 4 / 3 

was chosen to represent the Fowler-Hoyle model. The mass was 
ten times that of the sun and the initial radius was R = lO^cm. 

The Chiu model was approximated by an isothermal core 
of 1.435 Mq and a y “ 5/3 envelope giving a total mass of 
1.952 M 0 . This model might correspond to the centrally- 
condensed core of a more massive star with giant structure. 
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the extensive envelope being neglected. Although tuis is not 
as centrally condensed as the model with a positive temperature 
gradient in the core which Chiu first suggested, it does give 
a highly condensed structure and is in sufficient contrast 
with the Fowler-Hoyle model. 

B. Dynamical History of 10H -, Pre-supernova Models of 

Polvtropic Structure. 

In order to clarify the effects of neutrino energy 
transfer, the following models are to be presented: 

(1) the "no-neutrino" model in which all neutrino energy 
transfer is ignored, (2) "neutrino sink" model in which 
all neutrinos, once formed, are assumed to escape the star 
without interaction, (3) "neutrino diffusion" model in which 
diffusive energy transfer by neutrinos at high densities and 
temperatures may occur . The choice of a lOM^ polytrope of 
index 3 allowed comparison with the results of Colgate and 
White (1964) . The initial evolution of all three models 
(IOMq polytropes of index 3), was identical, following the 
path ABC in Figure (4) ; consequently this part of the 
evolution will be discussed only once. 

Figure (4) illustrates the history of one representative 
zone, (M r = 1.5^) falling into the core, of the models just 
discussed. The contraction was initiated by introducing a 



EVOLUTION OF REPRESENTATIVE 



LOG.J DENSITY- GM / CM 
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small energy sink to evolve the model. This was done by 
uniformly increasing the combination pair-annihilation and 
plasmon neutrino rates throughout the model until it con- 
tracted to the collapse point, slow enough to follow a 
Y - 4/3 adiabat more or less, but fast enough to do so 
without using too much computing time. The stability condi- 
tion (Appendix equation C. 1) on che hydrodynamic difference 
equations i' .nduly restrictive for the quasi-static 
evolution of hydrostatic models, so that this accelerated 
contraction was necessary. So long as the configuration 
follows the correct path in the prossuro-tompor.il ure plane 
for pre-collapse, and at the point of collapse has negligible 
kinetic energy, this approximation is valid. This part of the 
evolution carried the zone along the path AB. 

The evolution along the path BC represents the evolution 
of the material through the Fe-He phase change. At point B 
energy is removed by endoergic nuclear reactions so that th~ 
zones begin to fall inward rapidly. At point C the conversion 
of the iron-peak nuclei to helium and nucleons is nearly 
complete. In order to keep the amount of data manageable, 
quantities were printed out every 200 time cycles of the 
hydrodynamic calculation. The evolution was so rapid at this 
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point that the accurate position at which the transition was 
completed is not known, but it is probably as drawn in Figure (4). 

1. "No-neutrino" model. 

The evolution from C to D in Figure (4) may be explained 
as follows. The implosion proceeds until the nucleon terms 
in the equation of state provide sufficient pressure to halt 
the inf ailing material. When neutrino energy lo: i. ; : noglerlod, 
the following artificial situation develops. When the electron 
Fermi energy rises sufficiently to cause inverse beta-decay, 
the number density of electrons does not rise much with an 
increase in density, and hence the pressure contribution due 
to electrons does not rise either. The free nucleons which 
are formed do contribute an ideal gas term to the pressure 
(in this model). Since neutrino energy loss is neglected in 
this model, the material falling upon the initial 0.5M Q core 
is shock-heated to high temperatures. 

The path CD describes the thermodynamic history of the 
zone as it encounters this stationary core shock. The kinetic 
energy which the zone gained upon falling to this density 
was converted to thermal energy of an ideal nucleon gas. This 
i may be seen from the following estimates: 

I 

I 

I 

i 
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Change in potential energy upon 

contraction ^ -5.3 X 10 52 ergs 

Change in thermal energy necessary 
to support the core material 

hydrostatically ~ 4 . 2 xl 0 52 " 

Energy lost in photodisintegration 

of nuclei ^ l.OxlO ^ 2 " 

Energy left. to for in an over- ^2 

pressure and mass ejection £ O.lxlCT " 

The energy available to form an overpressure is negligible as 

far as an explosion of the model is concerned. In fact the 

core continued to adjust itself as overlying layers continued 

to rain down, and the implosion was not reversed. During this 

period the zone plotted in Figure ( 4 ) evolved along path DE. 

The calculation was terminated when a core of about three 

solar masses was formed; at this point there was no indication 

of any possibility of mass ejection. Although the histories 

of only "representative" zones are given in the Figures, 

v statements of results (as here with the absence of mass 

ejection) are based on examination of the behavior of all 

zones in the model. 

During a homologous contraction (or expansion) every mass 
zone of hydrostatic gas mass has its pressure related to its 
density by 

P = P (£_ ) 

0 p o 


4/3 
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where P q is the pressure and p the density of this zone for 
some reference configuration. This relationship defines the 
path in the pressure-density plane upon which the mass zone 
must lie for hydrostatic configurations (which are homologous 
to the reference configuration) . When the gravitational 
acceleration 

g(r) = - G V' r > 

r z 

is linear in the radius r, a gas sphere in free-fall contracts 
homologously. This occurs when 

M(r) - r 3 

or equivalently, when the density is constant. 

p = constant . 

For a density which decreases with radius, the outer zones 

of the gas sphere must be accelerated less than would be 

necessary for homologous contraction, and "left behind". 

The innermost zones of a gas cloud freely falling under its 

own gravity tend to fall homologously, leaving behind those 

4/3 

zones which do not. Thus the relation P -v p defines a 

locus of hydrostatic configurations for the imploding core 
(and for the whole star model in so far as the homology 
requirement is satisfied) . This curve in the P-p plane, and 
its corresponding curve in the p-T plane, define a sort of 
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stability boundary, such that if the pressure predicted by this 
relation is greater than the model actually has for a given 
density, ther- the core is unstable toward contraction and 
cannot be static. The dashed line in Figure (4) defines such 
a boundary. Note that the evolutionary path of the "no- 
neutrino" model does not cross this line, although at point 
E the zone shown is a part of a quasi-static core of 3M 0 
and the stability boundary is approached. By this time the 
structure of the core is by no means a scaled-down version of 
the structure of the same matter in the initial hydrostatic 
model. Although the homology requirement for the use of the 
stability boundary is violated, in fact there is no evidence 
that any significant reversal of implosion or mass ejection 
is likely. For smaller mass cores, it will be seen that the 
stability boundary concept is useful. 

2. "Neutrino Sink" Model. 

The "neutrino sink" model differs from the "no-neutrino" 
model just described through the inclusion of an energy loss 
rate of the form: 

Q = - 4.3 x 10 15 _ 1.1 (T 9 3 ) P erg/gm/sec 

P 

This is the sum of approximate forms for electron pair- 
annihilation and plasmon decay neutrino loss processes 
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mentioned earlier. It is not maintained that this expression 
is the correct form to use. Quite the contrary, other 
processes such as the URCA process for nucleons are probably 
more important. What is important is that the neutrino 
energy production rate is almost certainly as large as this. 

If it is assumed arbitrarily that all neutrinos, once formed, 
escape from the model, then this energy loss rate is a lower 
1 imit . 

What is the purpose of such an artificial model? Simply 
this: it shows that for this comparatively mild energy loss 

rate, there is no reversal of implosion, no mass ejection, 
but just the accumulation of a degenerate core of ever- 
increasing size when neutrino diffusion is neglected. This 
is shown in Fig re (4) , where this "neutrino sink" model 
follows the path CFG. In the segment CF the graphed zone 
encounters a stationary core shock, in which the infalling 
zones are slowed and become part of the core. For this model, 
the neutrino energy loss rate keeps the temperature much 
lewer than was the case for the "no-neutrino" model. This 
causes the core to form at a much higher density than was 
the case in the "no-neutrino" model. This means that more 
gravitational potential energy has been released, but is 
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lost from the star. The continued loss of energy Hows 
the core to evolve slowly to higher densities, along the 
path FG. The irregularities in the path along FG correspond 
to oscillations* cf the core, perturbed by the continued 
infall of matter . The mattor is quite degenerate at this 
point, so that temperature has little effect on the equation 
of state. Notice that the evolutionary path CFG never comes 
near the "stability boundary". 

3. "Neutrino Diffusion" Model. 

Neutrino opacity and the approximations involved in the 
assumption of diffusive energy transfer have been discussed. 

In this model the possibility of energy transfer by the 
diffusion of neutrinos is considered. It should be emphasized 
that the neutrinos were assumed to be in thermal equilibrium 
with the other particles in order to avoid a kinematic 
calculation. 

The temperatures shown in Figure (4) are unreasonably 
high because of the approximate nature of the equation of 
state. In particular, the thermal contribution of the 
nucleons to the pressure was underestimated. The neutrino 
opacity depends sensitively upon the average neutrino energy, 

* These oscillations may be due to the finite size of the 
mass zones as well as to the excitation of actual physical 
oscillations. 



38 


which depends upon the temperature. Because of this over - 
estimate of the temperature, the collapse of a model using 
the electron-neutrino scattering opacity behaved much like 
the previously discussed "no-neutrino" model, that is, a hot 
core was formed, but there was no mass ejection. In order 
to examine the effects of neutrino energy transfer, the 
opacity in the core was kept low enough so that the ^-.ergy 
transfer time scale was of the order of the hydrodynamic 
time scale. This affects only the zones of hiq ar density . 

p ~ 10 12 gm/c" 3 

so that opacity in the crucial region in which the infall 
of the matter is reversed, is just that for noutrino-olcctron 
scattering, with average neutrino energy E ~ 3kT. A report 
on investigation of the validity of this assumption for the 
opacity at high densities is currently being prepared for 
publication. 

The "neutrino diffusion" model was calculated in 
exactly the same manner as the "neutrino sink" model until 
the neutrino mean free path became short compared with the 
dimension of the star. At this point the transfer of energy 
by neutrino diffusion was calculated, and some results are 
indicated in Figure (4) . The path CD is identical with the 



"no-neutrino" case, but the temperature continues to rise 
rapidly as DH shows. 

That this chock heating i:; more pronounced f.hon in the 
"no-neutrino" model is quite important. Why this is the 
case may be seen as follows. Neutrino diffusion removes 
energy from the core. Before the neutrino diffusion is 
initiated, neutrinos are unhindered as they escape the core, 
again removing energy. This loss of energy prevents the 
temperature from rising rapidly, so that it is nucleon 
degeneracy pressure rather than thermal pressure which halts 
the infall of the core. This means that the core will have 
a much greater density than in the case of the "no-neutrino" 
model, which in turn implies that more potential energy 
’ s released by the contraction. Hence more energy is 
available for expelling matter. 

The situation is now unstable in the following sense. 

If the infalling matter supplies kinetic energy to the core 
faster than this energy can be removed by neutrino diff sion, 
the temperature will rise. Because the neutrino interaction 
cross-sections are roughly proportional to the square of 
the neutrino energy, and because at higher temperatures 
neutrinos are formed with higher energies, the opacity 
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increases. The transfer of energy by neutrinos then decreases, 
and the temperature rises still more. Thus the medium may 
become opaque if the inflow of kinetic energy is sufficiently 
high. 

The greater potential energy released in the "neutrino 

; 

diffusion" model is now available for reversing the implosion 
of the outer layers. A neutrino diffusion wave sweeps 
out of the core, leaving tie matter behind it opaque. The 
zone shown in Figure (4) is heated so that the path DH lies 
well above the stability boundary. Along the path HKL, this 
zone falls into the core. The path KL shows the core adjusting 
itself while overlying matter is ejected. Figure (5) shows 
the curve already shown on Figure (4) as CDHKL but is now 
labeled MNOP. The new curve on Figure (5) , DBF, corresponds 
to the zone with an interior mass 

M r e 2.0 Mq 

and which is ejected from the star. The path DE corresponds 
to the heating of the zone by the neutrino diffusion shock 
wave, and EF is the subsequent expansion of the zone as it 
leaves the star. The remnant core mass was 1.8 Mg. 

C. Comparison With the Results of Colgate and White. 

Although the hydrodynamics of supernova envelopes has 
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received much attention, the only previously published 
investigations (of which the author is aware) of the initial 
instability and subsequent dynamical history of supernova 
interiors are those of Colgate and his collaborators. A 
review of these investigations by Colgate and White (1964) 
describes some of their more recent rest Its and is the primary 
source for the following description of their work. In 
particular, the evolution of type II supernova models was 
followed by means of a numerical hydrodynamic computer code 
from the onset of gravitational collapse to the reversal 
of the infall of the core (due to terms in the equation of 
state corresponding to a nucleon hard-core potential) . The 
loss of neutrinos emitted in inverse beta-decays cools 
the core during implosion. This loss rate is approximated 
by a simple analytic form; a more exact analysis would 
involve the evaluation of Fermi-Dirac integrals because the 
material becomes degenerate. A partial deposition of this 
neutrino flux in the stellar envelope and the shock wave 
reflected upon the formation of a neutron star core provide 
sufficient energy to eject ~ 80% of the mass of a 10 solar 

mass star. To simulate the emission and deposition of neutrinos 

\ 

from the shock wave formed by the infalling material raining 
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down upon the quasi-static core, half of the emitted energy 
was deposited in the matter external to the core shock. 

This deposition was initiated only when the core shock was 
formed and was turned off when the rarefaction due to ex- 
pansion terminated the core shock. The time-dependent 
energy sink term, integrated over the core, is 



where the factor in parenthesis is just an analytic 
approximation to the inverse beta-decay neutrino loss rate, i.e. 

d g y * - 0.1 T erg/gm/sec 

dt 

3 

where p is the density in grams/cm and T the temperature 
in KeV. The rate of neutrino energy deposition that was 
used, in units of ergs/gm/sec, is 

a -*■ a m exp (-K f r p dr) 

W S deposited 4nr 2 xp Jr^ 

for 

r r 

shock 

and where 

0 » 

K = In 2/(1* p dr 
.shock 


r 


m i m m 


n i 
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The test for initial equilibrium proceeded for a real time of 
30 seconds for a 10M o polytrope of index 3 taken to be the 
model of a pre-supernova star. This model was inspired by 
considerations of Fowler and Hoyle (1964) . The instability 
was initiated by removing 1% of the internal energy. The 
core formed cold with 5% of the mr»ss of the star; after the 
implosion was reversed in the innermost mass zones a shock 
formed and neutrino deposition of energy was initiated. In 
this case 2M© accumulated in the core before sufficient 
heat was deposited to reverse the implosion of the outer 
layers and create an explosion. 

Colgate and White (L964) also discuss some calculations 

involving initial models of 2.0 and 1.5 solar masses. These 

models evolved on such a low adiabat that, rather than pass 

through Fe-He phase transition, they were brought to dynamic 

implosion by rapid electron capture (inverse beta-decay) at 

1 3 

a density above 2x10 “ gm/cm . The subsequent core formation, 
shock wave, neutrino emission and deposition, and finally 
explosion followed as in the 1OM 0 case. The expansion 
velocities and residual core mass were lower, but not 
drastically so. 

Figure (6) compares the result of the "neutrino diffusion" 
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model to the results of Colgate and White (1964) . Although 
the data for the plot of Colgate's model were taken somewhat 
crudely from graphs, there seem to be three differences. 

This first concerns the Fe-He phase change region BCD which 
lies on a lower curve for the diffusion model. This does 
not affect the dynamics much since both models are in free- 
fall at this point. Considering the widely different methods 
used in treating this phenomenon, the difference in the two 
paths is not surprising. The second point is that Colgate's 
sudden heating (path DE and MN in the diffusion model) seems 
to occur at lower densities. This is thought to be attributable 
to Colgate's technique of depositing energy. The last difference 
is the rapid cooling of the core as shown by Colgate's model. 
Colgate's energy transfer technique probably is inaccurate 
at this point; this cooling occurs after the mass ejection 
so that its effect upon other aspects of the supernova 
phenomena is small. Also, the diffusion approximation will 
incorrectly predict the energy loss rate as the distribution 
function for neutrinos departs from its form for thermal 
equilibrium. The problem of core cooling would be properly 
handled only by a detailed transfer calculation. 
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D. Effect of Pre-collapse Structure, 

Figure (7) illustrates the history of two roproscntotivo 
zones of the centrally-condensed pre-supernova model described 
previously. This model was evolved slowly along a path 



until electron capture instigated collapse. This phase, 
similar to the evolution of the 101^ polytrope model before 
the Fe-He phase transition, Is shown in Figure (7) as the 
paths AB and MN. By the time the ? nones had reached the 
vicinity of points B and N, the ccie had "bounced" and the 
neutrino diffusion wave was moving outward. 

The effect of this diffusion "wave" may be seen in 
Figure (7) . Consider the zone with 0.98 Mq underneath it 
first. Along path NP there is a wiggle which was caused 
by initiating the diffusion calculation and has no significance. 
During this time the diffusion wave has not yet reached the 
zone. Path PQ shows the heating of this zone as it, falling 
in, encounters the diffusion wave sweeping outward. During 
the time the zone moves from Q to R it is inside the 
neutrino emitting surface. It falls on the core at point R, 
and evolves slowly from this point on. Further evolution is 
due to "slow" energy loss by neutrino diffusion (the time 
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-3 

involved between R and T is t ^ 10 sec) . 

The ejected zone = 1.48 f^) encounters a combina- 

tion of shock heating and neutrino diffusive heating along 
the path BC. At C, the peak temperature for this zone, 
expansion begins. Along CD there is some heating due to 
acceleration of lower zones, but along path DE the pair- 

annihilation neutrino energy losses make the path slightly 

3 

steeper than p ~ T . At point E, the thermally decomposed 
nuclei begin to recombine by exoergic reactions. This 
causes the temperature to drop off more slowly along path 
EF. At this point the calculation was terminated. 

The similarities and differences in the two structur- 
ally different models are summarized in Table I. The 
"velocity of ejected matter" quoted in Table I means the 
average velocity of that matter behind the ejection shock 
wave at the conclusion of the calculation. The velocity 
corresponding to the observed expansion velocities of 
supernova remnants should be less because an envelope of 
several solar masses still lies outside the ejection shock. 
Extension of these calculations to this asymptotic ejection . 
velocity is contemplated. 

The more violent ejection of matter from the centrally 


condensed model is to some extent due to nuclear recombi- 
* 

nation late in the expansion. The similarities are more 
striking than the differences. These two widely differing 
models, brought to collapse by different mechanisms, never- 
theless have surprisingly similar characteristics. The 
remnant core mass seems to be higher for the le ,s condensed 
model; this might be expected because the mass of the sphere 
inside which the gravitational acceleration. 


g(r) 


GM(r)_ 


is nearly linear, is smaller in the more condensed model. 
This determines the mass of the material which halts its 
contraction as a unit; this material collapses approximately 
a uniform density sphere would. 

E. Comparison of Calculations and Observations . 

Because of the low frequency of occurrence of supernova 
outbursts (about 1/century /galaxy) observational information 
is meager. Zwicky (1956) gives a history of supernova 
observations. More recent accounts by Shklovski (ISoO) 

and Minkowski (1964) make it reasonable to identify type II 

\ 

supernovae with the catastrophic disintegration of a massive 


* The polytrope model was not followed this far. 
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star Table IT contains a summary of some character- 

istic properties of supernovae type II. Details concerning 
the optical spectra and light curves, which are not investigated 
here, may be found in the references mentioned above. 

Comparing Tables I and II, it will be seen that there is 
reasonable agreement between the observational evidence and 
theo' 'tical predictions. The visual magnitude expe' -ad from 
the theoretical models has not been estimated, but the kinetic 
energy of ejected mass is appropriate. The velocity of ejec- 
tion appears larger than observed, but in fact the asymptotic 
value of the velocity will be lower, especially if the star has 
an extensive envelope as has been assumed. This comment applies 
to the condensed model, in particular, which is envisaged as 
having a large envelope. 

IV, Interpretation of the Results of the Model Calculations . 

A. Peak Temperature and Muon Neutrino Energy Los s. Muon 
neutrinos do not have the same interaction cross-sections as 
electron neutrinos. Unfortunately the opacity for mu neutrinos 
in a hot, dense medium like a supernova core is not as well in- 
vestigated as for electron neutrinos. It appears that neutrino 
production by muon-pair annihilation will dominate production 
from simple muon decay, for higher temperatures at least, as will 

be shown by the following argument. 

# 

Some rough estimates of the muon neutrino energy loss rates 
may be made as follows. Ifre number density of fermiou pairs. 


neglecting any decay modes , is given by 


2 

“ ■ h 3 r tt 


d V 


exp [ (C + +ji) AT] 


Assume that the muons are nondegenerate and that the number 

* 

densities of |i • and u are equal. Then, since 


exp C (C +fi ) /kt »1 

where the maximum of the integrand occurs, then 

. * p ■ Or) 3 V f (8) 


2 . 


where is the number density of muon pairs, me is the 
rest energy of the muon, and 


CD 

f(p) = f esqp C-b (1+x 2 ) 1/ ^ 2 ] x 2 dx 

o 

2 


B = 


For large 3, that is, kT «mc , 


me 

kT 


-0 


f(e)«/Tr/2 


so that 


IQ 38 

N p ~ 3.2 


,V2 


The rate of neutrino energy loss at high enough temperatures 
by muon decay is roughly 


& ~ - 2Np C/(pT)erg/gm/8ec 


* Chiu and Stabler (1961) present this approximation for 
electrons . 
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where 5 is some average energy of the emitted neutrino, p 
the matter density, and t the half life of the muon. 

Assuming equipartition of energy, £ might be of the order 
of 35 MeV. 

Using the same approximations of nondegeneracy and 
2 

kT « me , Chiu and Stabler give the approximate pair 
annihilation energv loss rate 

3 -(2m c 2 AT) 

18 a q e e 

s e » - 4.8 x 10 ^ ergs/gm/sec 

Correcting this for the heavier mass muon, 



- (2400/T g ) 


Neglecting the effect of muon decay on the number of 
muon pairs formed in equilibrium with radiation may be 

justified as follows. The characteristic time for muon 

-6 2 
deacy is 2 x 10 seconds. Taking kT « m^c the muon pair 

annihilation cross-section has roughly the same magnitude as 




2 

cm 
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So that the mean reaction time is 


T ~ 


r N a v 


which for T = 120 gives 
9 


, -13 

t ~ 10 sec 
r 


so that 


t « t decay 

and this implies that muon decay does not greatly alter 

the muon number density. Some estimates for the energy 

loss rates by muon type neu . * tos are given in Table III. 

The energy loss rates have units erg/gm/sec and are 

12 3 

all evaluated at a density of 10 gm/cm . The relaxa- 
tion times for these rates, if they proceed unimpeded 
are given in Table IV. 
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In the last case the energy Loss by neutrinos emitted 
by muon decay is dominant, but the time scale for cooling 
is then longer than the characteristic collapse and mass 
ejection time scale. From these rough estimates it might 
be expected that the temperatures for the n=3 polytrope 
model. Figure (5), aie artificially high. The muon- 
neutrino energy loss is probably not so large for the 
isothermal model. It is not convincing to estimate what 
effect muon neutrinos will play in supernova explosions 
without a careful analysis of all possible reactions, 
both for neutrino production and neutrino opacity, and a 
careful estimate of the temperature. It does appear that 
larger remnant cores might be expected with the inclusion 
of muon- type neuter.' no energy loss, but little more can be 
said at this time. 

B. Electron Neutrino Luminosity and Detectability . 

The immense energy radiated by neutrinos that the previous 

models preset, and the high temperature of the emission 

surface, suggest that it might be possible to detect super- 

* 

novas by their neutrino flux. Dr. Raymond Davis, Jr. of 


* 


R. Davis, Jr. (1965) 
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Brookhaven haa kin-oily provided information about his detector 
with which to evaluate this possibility. 

If a detector contains N absorbers, then the fraction of 
the total integrated neutrino flux absorbed in the detector 
is 



where c is the interaction cross-section and R is the dis- 
tance of the source. The total integrated flux is roughly 

E 

F . -**• 

* 3kT 


where E tQt is the total energy emitted in the form of 
neutrinos (only electron neutrinos are considered because 
of the uncertainty in the muon reaction rates), and T is the 
temperature of the zone just inside the neutrino emission 
surface. 

For kT ~ 10 MeV, and 


i, in 53 

E tot. ~ 10 ergs 
N ~ 2 x 10 30 Cl 37 atoms 

o ~ 7 x 10“ 41 cm^. 


where the cross-section is for the reaction 

,37 


V + Cl' 


♦ A 37 . 
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and for a minimum detectable signal of .100 counts, 

22 

R ~ 1.5 x 10 cm 
~ 5 kpc. 

the maximum source distance is 5 kiloparsecs. According 

* 

to Allen , the diameter of the Galaxy is 

D ~ 25 kpc. 

Recalling the supernova rate quoted earlier, one per century 
per galaxy, it appears that no burst of electron neutrinos 
is likely to have produced a detectable signal since suit- 
able detectors have been operating. 

C. General Relativity and Core Collapse . 

In order to check the neglect of general relativity, 
the Schwarzschild radius 


2 GM 


8 


0.3 x 10 


6 c — ^ 


.cm. 


c 

is compared to the radius of the dense core. 

Condensed Model 


n = 3 polytrope model 


r /r core 0.22 0.38 

s 

As anticipated the general relativistic effects are becoming 
important, especially for the n * 3 polytrope model. These 
effects will be more pronounced in the core. The mass 
ejection, occurring before such high densities are reached 


* Allen (1963), p. 267 
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V# 

(r /r m* "r>z the n = 3 polytrope is then only ~ 0.16), will 

S CO T ; 

be less sensitive to this effect. However, the bounce point 
for the core collapse, where the infall is halted, is some- 
what dependent on both the nuclear equation of state and the 

temperature of the infalling matter. ' Careful investigation 

/ 

with general relativistic hydrodynamic equations is needed 
to confirm that neglect of general relativity is reasonable 
until after mass ejection. 

V. Critique of Calculational Method . 

Approximating the transfer of energy by neutrinos by a 
model of diffusing neutrinos in local thermodynamic equilib- 
rium is probably correct in a rough sense, but while it may 
be on firmer theoretical ground than the intuitive approach 
of Colgate and White (1964) , it leaves much to be desired. 

On the other hand, the solution of a transport equation 
coupled with hydrodynamic motion poses extremely difficult 
computational problems. 

The shock width in Figure (5), path MN, and in Figure 
(7), path PQ, stretches over a density range of more than 
a power of ten. Direct examination of the numerical 
results of the calculation reveal that in general the 
width of the shock zone is smeared over too wide a 
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region due to the lack of an adequately fine zoning mesh. 
Unfortunately, the present generation of computers is too 
slow to make the use of models with more zones feasible. 

The computing time goes as the square of the number of 
mass zones in the model, so that this limitation is diffi- 
cult to overcome. 

By the similarity in energy release and core size as 
obtained by the polytrope of index 3 and the isothermal 
model, it appears that the interior dynamics of supernova 
are relatively insensitive to the structure of the pre-super- 
nova model. This does not mean that the existence of an 
extensive envelope such as found in massive red giants might 
not affect the velocity of ejection and the peak shock 
temperature in the matter ejected. It does mean that the 
implosion and "bounce '• of the core, as well as the neutrino 
energy transfer process are insensitive to initial structure. 
The interior dynamics is much the same for the two models 
presented. 

Summary 

The calculations reported here indicate that it is 
possible to construct reasonable models of supernova by 
assuming energy to transferred by electron- type neutrinos 
during stellar collapse. Whether or not considerable mass 


/■ 



ejection by this process actually occurs depends critically 

upon the average neutrino energy and the opacity for neutrinos 

12 3 

in regions of high density (p > 10 gm/cm ). Unfortunately 
these quantities remain uncertain. Because of the temperature 
uncertainty, the emission rate of muon- type neutrinos is also 
unknown. Even the relatively low estimates of energy loss 
rates made for muon-decay (neutrinos from pion-decay may be 
more important) indicate that muon- type neutrino processes 
will be of paramount importance at the higher temperatures. . 
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Appendix on Numerical Methods 


Because the techniques used in this research are not 
as yet well known to most physicists and astronomers a 
brief summary is presented here. The equations of hydro- 
dynamics may be written as follows: 

1) mass conservation 


dM R = 4TTR 2 pdR # 

2) momentum conservation 



3) energy conservation 

dE - s dt - PdV, 


but using 


gives 


dE 





(« at - (p + (|f ) ) av ), 

T 


4) equation of state as a function of temperature T and 
specific volume V 


P * P(v,T) , 

m ET(V,T), 


r 
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EV (V, T) . 


T 

In these equations is the mass interior to some radius 
R, p = 1/V is the density, U is the velocity of Lagrangian 
mass element at R, E is the internal energy per unit mass, 
s the rate of addition of energy per unit mass, and P is 
the pressure. 

A. Difference Equations . As they stand, the fluid 
dynamic equations are highly nonlinear. Because of this 
difficulty only few analytic solutions are available, 
and all of rather limited applicability. For ether cases 


one is usually forced to some sort of approximation tech- 
nique which is often as complex as numerical solution, 
and which may tend to obscure the physical situation. 

Even the approximation techniques are generally restirc- 
tive. In view of these problems it is often the case that 
numerical solution difference equations is preferable. At 
first we neglect radiative transfer of energy. Neutrino 
energy transfer will be treated subsequently. 

The star will be divided into concentric spherical 
shells by J boundaries numbered 1, 2, 3,...J (from the center 


outward) . Quantities associated with the zone boundaries 
will be subscripted j ; those associated with zone centers 
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are subscripted j + 1/2. Time centering is indicated by a 
superscript n in a like manner. 

After using some simpler forms, difference equations 
quite similar to those of Colgate and White (1964) were 
adopted. Although they have been presented in the above 
reference, the equations are discussed here for complete- 
ness. 

The initial configuration is input. The position 
and velocity of each boundary 

R 1 
D 


j = 1, • . . JJ 


ir 


where JJ is the number of boundaries, the pressure, specific 
volume and temperature of each zone 

P j + 1/2 
V j + 1/2 

T 1 

* j + 1/2 

must be specified. The mass of each zone may be calculated 
from 

.1 


DM 


j+1/2 


f (*j+0 - Of! ** 


j.' ' j+1/2 


A.l 


b 

i 
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The effective mass of a zone as seen from a boundary is 
the simple average 

DH j = * C m Ul/2 + DM j-l/2) 

A. 2 


which assumes that the zone masses differ little. In 
practice it was found that results were often better when 
neighboring zone masses were the same or changed by a 
small constant fraction. The total mass inside a boundary 
j+1 is 

XM 1 , = XM 1 + DM'!' . 

D+l j 3+1/2 


A. 3 


In the Lagrangian system mass is necessarily conserved 
until the zoning is changed. When a system changes drasti- 
cally it is sometimes possible to shorten the time required 
for a calculation by rezoning the configuration, but other- 
wise the zone masses remain the same. 

The equation for momentum conservation may be written 


as 


yn+1/2 _ y n- 1/2 

j 3 

A. 4 


,„n.2 rr n „n , _n-i/2 n-1/2"] A n 

- ( V L^j +1/2 - P j-l/2 + Q j+l/2 • Q j-l/2J 4t /DM j 


G 


2l 


( r ") 2 


At 


n 
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where the fluid velocity is 


U 


dR 

dt 


so 


A. 5 


„n+l „n n+1/2 .,n+l/2 

R. = R. + U. ' At ' 

3 3 3 


From this the specific volume car be updated by 


v n+1 7 

3+1/2 


2 

3 


< R j + ^> 


n+1. 

~ (R j } 


“* 3 + 1/2 


A. 6 


which reflects mass conservation. The specific volume 
evaluated at the same point in time as the fluid velocity 
U, that is at n+1/2, will be useful. 


♦ ^+i /2 > 


A. 7 


At this point a linear extrapolation in time is made for 
the new temperature at point n+1/2 . Initially 


3/2 

j+1/2 


but afterward 


T n +-/2 _ n 
3+1/2 3+1/2 


1 

+ 2 


At 

At 


n+1/2 

n-1/2 


A. 8 


T 


1 

j + 1/2 


(T 


n 

j+1/2 


i n " 

3+1/2 


) 


This will be used to determine the temperature at the epoch 
n+1. 
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Evaluating the equation of state at time step n+1/2 


gives 


A. 9 


n+1/2 . 

P j+l/2 = P(T 


n+1/2 n+1/2. 

j + 1/2 ' J + 1/2 * ' 


A. 10 


A. 11 



n+1/2 

) 

V j+1/2 


ET (T 


n+1/2 
j + 1/2' 


vn+i/2) 
J + 1/2 




n+1/2 

) 

T j+1/2 


= EV 


, n+1/2 
U j+1/2' 


vTK 2 2 ) ■ 

jfx/2 


The quantity Q in the momentum conservation difference 
equation is the so-called pseudo-viscosity term which 
stabilizes this set of difference equations. When zone 
boundaries approach rapidly it supplies a large pressure 
to prevent them from crossing. In a shock the pseudo- vis- 
cosity term converts kinetic energy of zone motion into 
thermal energy, and is negligible elsewhere. The form 
used is 

n+1/2 2 ftjN+1/2 _ if 

j+1/2 ( J+l U j 'j+1/2 lf J+1/2 J + 1/2 


. ..n+1/2 
and u j+ i 


< u n+1/2 
3 


otherwise, 


. A. 12 
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which is zero on expansion. 

The energy conservation equation becomes 

T n+1 = T n + f _ (P n+1 / 2 + 0 n+l/2 n+1/2 

T j+l/2 j + 1/2 + n+1/2 L 1 j+1/2 + Q j+l/2 J+1/2' 

et :+i/2 


A. 13 


(v n+1 , 

V j+1/2 


^+ 1/2 } 


+ 


n+1/2 A n+1/2 
^j+l/2 4t 


Although the energy source term s has not been specified, a 
form for it could have been evaluated at epoch n+1/2 along 
with (A. 9) for instance. 

B. Pseudo- viscosity Technique . The pseudo-viscosity tech- 
nique for treating hydrodynamic shocks is due to Von Neumann and 

* 

Richtmyer (1950). There are few references to it in the litera- 
ture although it seems to be arousing some interest among 
astrophysicists. 

Attempts to solve the equations of fluid motion by numeri- 
cal procedures are greatly complicated by the presence of shocks. 
Mathematically the shocks are represented by surfaces upon which 
the temperature, density, pressure and fluid velocity are dis- 
continous. The partial differential equations governing the 
motion require boundary .conditions connecting the values of 


* 


The author found Richtmyer (1957), Fromm (1961), Henyey 
(1959), Christy (1964) and Colgate and White (1964) most 
useful. 
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these quantities on eaich side of the shock surface. The 
Rankine-Hugoniot relations, i.e., local conservation of 
mass, momentum and energy by the fluid, supply the 
necessary restrictions, but are difficult to apply 
during a calculation because the shock surface moves 
relative to the mesh points in space- time which are 

\ 

used for the numerical work. The nonlinearity of both 
the differential equations and the boundary conditions 
does not simplify the problem. The motion of the shock 
surfaces is not known in advance but is determined by 
the differential equations and the boundary conditions 
themselves . 

The method of Von Neumann and Richtmyer automatically 
treats shocks and avoids the necessity for pre-knov»ledge of 
shock motion by utilizing the effects dissipative mechanisms 
(such as radiation, viscosity, and heat condiction) upon shocks. 
When viscosity is considered, the mathematical shock discon- 
tinuity becomes a thin layer in which the pressure, density, 
fluid velocity and temperature vary rapidly but continuously. 

By introducing an artificial dissipative mechanism to spread 
this shock layer over a few mesh points, the difference 
equations approximating the equations of fluid motion can be 
used throughout the calculation, as if no shocks were present. 



7 

? 
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In the numerical results the shocks appear as rapid changes 
in the variables, almost discontinuities, which have very 
nearly the correct speed and across which the pressure, 
temperature, and density have very nearly the correct changes. 

In an actual physical problem the dissipative mechanisms 
are generally much smaller than the artificially introduced 
viscosity term. The limit on computational reproduction of 
a physical situation is that the zone size be smaller than 
the smallest dimension of interest. The quadratic depend- 
ence of (Q) on the velocity difference insures that this 
form for the artificial viscosity is small except in the 
shock region. Note that the 1/V dependence gives an in- 
creasing pseudo- viscosity for large compression. 

C. Stability of the Difference Equations . For a more 
complete discussion of the stability of finite difference 
approximations the reader is referred to Von Neumann and 
Richtmyer (1950), Fromm (1961) and Richtmyer (1957). To 
clarify the meaning of stability, consider the exact solu- 
tion Y(r,t) to the one- dimensional differential equations 
of fluid dynamics for some specified initial- value problem. 

Let ¥]? be the corresponding solution to a system of differ- 
ence equations which approximate these differential equations 
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The error of the approximation is then 
j Y^ - Y(r=jAr, t=n^t) | 

Suppose that this error is small enough at some time t, and 
consider some small perturbation 6 Y . Does this perturbation 
grow with time? If so the difference scheme is obviously 
unacceptable as an approximation to Y(r,t). 

The criterion for stability against such small per- 
turbations, for a set of difference equations such as pre- 


sented here, is that the time step for integration At satisfy 

At 


C.l 


< * x i+i/2 

v 

s 


where Ax is a zone w idth and v g is the local sound velocity. 
This must be true for each zone. There is a simple physical 
interpretation for this restriction. Neglecting any sort of 
radiative transfer and considering only gasdynamic effects, 
the minimum time for material at zone boundary j to communi- 
cate with material at zone boundary j+1 is just the sound 
traversal time t 

s 


t 

s 



AX i+V2 

V 

s 


At < t 

s 


The requirement 




is simply that conditions at boundary j + 1, say, at time 


epoch n+1 be physically independent of what transpires 
at boundary j at the same epoch n+1. Most physicists 
are more familiar with the analagous situation in relativity 
where the velocity of light plays the part here taken by the 
sound velocity. Although such a restriction is mathemati- 
cally required of difference equations such as these, this 
restriction is physically necessary only when gasdynamic 
motions having velocities of the order or greater than the 
sound velocity are encounteed. If this is not the case a 

different set of difference equations might be developed 

* 

which had a less stringent requirement on the integration 
time step. 

For a complex problem in which there are drastic changes 
from the initial configuration, considerable computational 
time may be saved by choosing the integration time step At 
as the maximum value consistent with the stability require- 
ment (C.l). In the supernova problem it was discovered that 
this requirement, while necessary, was not sufficient to 
reproduce the physical situation faithfully. In the adia- 
batic contraction of a gravitating uniform sphere, it was 
found that the analytic solution was not successfully 


* See Henyey (1959) for example. 
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? 


approximated unless some restriction like the following 
was used: 


AV" +1/2 « 0.02 v”. 

so that the fractional change in specific volume was only 
a few percent. Colgate and White (1964) use the following 
simple and effective restriction on the time interval. 

* .^n+1/2 


n 


0 02 * V 

At n+3/2 ^ V j + l/2 


n ,31-1 | 

j+1/2 - j^l/2 


C. 2 


In the difference for.u for the energy equation, an 
energy source term s appears. This allows energy to be 
added or removed locally, although as vritten not 

explicitly account for _.»ergy transfer between z. ;es. 

When energy is suddenly added or lost instabil... '. ,.es often 
result. Now the energy deusri ty E of a flul-^ -;y be written 
as 

E ^ PV 

V - 1 

where P is the pressure, V the specific volume, and v the 
"ratio of specific heats" which is constant for an ideal 
gas. Generally y is a slowly varying function. With this 
in mind, the following time interval restriction was used 


r i ; 

i i 


% 
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n 


At 


n+3/2 


0 . 02 P7 . At 

1 + 1/2 j + 1/2 

P n v 1 " p n_1 

j+1/2 j + 1/2 j + 1/2 


n+1/2 




C. 3 


ana found to be adequate. 

Using these three restrictions on • time step, the 

new time interval can be calculated. The sound velocity 
was found to be adequately represented by the simple 
approximation 

v * J2PV 
s 

3/2 i/2 n 

The time intervals At and At are input. Then ^t , 

which is .eeded in (A . 4 ) , the momentum conservation equation. 


is 


..n 1 ,.^ 1 - 1/2 

At = — (At + 




C. 4 


This procedure allows a small, conservative estimate of the 
time interval to be input which insures that the stability 
requir jments are not violated. The scheme than chooses the 
optimum time interval and rapidly approaches it. It is 
noted in passing that the reason for centering these differ- 

•k 

ence equations in time is to allow this calculated time 
step scheme to be used accurately with varying time steps. 


* Only the pseudo- viscosity Q in (a. 13) is not centered in 
time, a condition mitigated by the fact that Q is not a 
physical but a computational quantity, and that computa- 
tions using this Q are accurate. 
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D. Boundary Conditions . In the supernova problem 
the interior boundary was taken to be at the origin, so 
by symmetry 


r" - 0 


uj = 0 

for all time. This was not necessary; the inner boundary 
could have been at some distance R from the center of the 
star, acting as a spherical piston with velocity 

U* = f(t n ) 

* 

which in general varies in time. There are difficulties 
associated with zones near a piston in this sort of scheme, 
however, and care should be taken. 

Having specified the center of the star as the inner 
boundary, symmetry assures that pressure, temperature and 
specific volume are continuous through the origin and no 
inner boundary condition need be specified. The outer 
boundary is free to move, however, 

s 

in general, so that its motion will be determined by the pres- 
sure (and artificial viscous pressure) at JJ+1/2. These are 
• not calculated and must be imposed. If 

t 

* Christy, R. (1964) . 
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_A » _* * 

jJ+1/2 = ” J J-l/2 

n-1/2 = _ qH- 1/2 

U JJ+l/2 -JJ-1/2 

D.l 

then the total calculational pressure (P+Q) will be zero 
at JJ. This boundary condition was used. In this particu- 
lar problem the motion of the inner regions was most inter- 
esting so that the choice of surface boundary condition did 
not happen to be critical. 

E. Analytic Checks of Numerical Results . In order 
to test the validity of the numerical techniques employed, 
several problems for which exact analytic solutions exist 
were calculated. Some of these results are to be found in 
Colgate and White (1964), and they are reproduced in the 
author's thesis both as an argument for the validity of 
these particular difference equations and because they 
provide insight into the technique. The problems for which 
checks were made are: (1) a strong, plane shock propagat- 

ing through an ideal gas with Y = 5/3 and density decreasing 
with the -7/4 power of the distance, (2) a strong spherical 
blast wave in ah ideal gas with a density 



/ 
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(3) the adiabatic collapse of an ideal gas sphere of uniform 

density, and (4) the hydrodynamic motion (or lack of it) of 

* 

a gravitating gas sphere in hydrostatic equilibrium. There 
is excellent agreement between the numerical and analytical 
solutions in all four cases. 

F. Radiative Energy Transfer . If the diffusion 
approximation is valid the change in energy density due to • 
sources and time- dependence can only affect lengths large 
compared to the mean free path. Then, for conservation of 
energy. 


ff +V yE = s- Av-? 

where s is the energy generation rate (per unit mass) due 

to sources, E the energy density per unit mass, v the velocity 

f the source, p the mass density, and I? the energy flux. If 

source motion can be neglected, and macroscopic changes occur 

on a time scale much larger than the mean free time for the 

diffusing particles, then 

«. ,^c dE Am . 

Ps * " 7 (_ 3 dT 4T) 

Assuming sperical symmetry, this may be rewritten as 


dT L (r) 


dR 


4ttR 


where r is the radial coordinate and 


l 


* The exact solutions for (1), (2) and (3) are to be found 
in Burgers (1949) , Sedov (1959) , and Colgate and White 
(1964), respectively. r 
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L(R) 


r 2 

J 4rrR s p dR 


D = 


Zc dE 
3 dT 


which are the standard forms for radiative diffusion usee 
in quasi-static stellar models. 

If macroscopic changes occur fast enough so that work 
done by pressure forces PdV must be included in the energy 
conservation equation, then 


dE = s dt (R 2 j* - PdV 

R p 

where spherical symmetry is assumed, V is the specific volume 

* 2 

and a is the radiation constant. Using dM = 4npR dR, 


where 


F.l 

and 


F.2 


dE 


(s 


dL 

dM 


) dt - PdV 


L * 


2 2 
(4ttP ) 


ac 

3k 


d(T 4 ) 

dM 



is the Rosseland mean opacity. Rewriting the energy con- 
servation equation in terms of the temperature T gives 


F.3 


dT 


o 


- - 


(P + 77, 


dE 

av 


) dV 


■]/<“ 


&T 


* For electron-type neutrinos and antineutrinos in thermal 
equilibrium and in equal abundance, the radiation constant . 
is a (neutrinos) *=> 7 a (photons). 

8 
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where (~)dt is the energy lost by diffusive transfer, PdV 
is the work extracted by macroscopic motion, and s is the 
energy gain by other mechanisms. 

G. Difference Equations and Boundary Conditions For 
Radiative Transfer , Equations (F.l), (F.2) and (F.3) may 

be incorporated into the hydrodynamic difference equations 
discussed in the previous section. Comparing (F.3) with 
our earlier energy conservation equation suggests the 
following difference equation, by analogy with (A. 13), 


« n+ ^ — rp n . — ,.,A~ — j” _ /p n tl/2 n+1/2 

3+1/2 3 + 1/2 + prp n+l/2 L ' 3 + 1/2 + j+1/2 

ET j+l/2 


n+1/2 _ '^i+1 - 

(8 j+ l/2 DM j+1/2 


-L 


) 


G.l 


But then (F.l) becomes 


.. n+1/2 l6TT^ac ,_n+l/2» 4 
“j — (Rj ) 

G. 2 


. n+1/2. 4 _ n+1/2 4 

( T i+l/2 j ir i-l/2 j 

(DM * AK)” +1/ ^ 2 


where 


(DM * AK) 


n+1/2 « I /mi __.n+l/2 AK n+1/2 ^ 

j 2 (DM j+l/2 ^+1/2 + m j-l/2 ^-1/2' 



.G. 3 



I 
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and (P.2) becomes 


AK n+1/2 = AK <t " +1/2 v" +1/2 i 

3+1/2 “ (T j+l/2' Vl/2' 


The term AL^ + ^ ' involves quantities evaluated at space points 

j-1/2, j. and j+1/2. The latter will not have been evaluated 

when T j + i/2 to calcu l at ©d from (G.l), if the method of 

sweeping through the space-time mesh described previously is 

used. The difficulty may be avoided by evaluating 

TT n +l/2 _n+l jh- 1 ..n+1/2 n+1/2 _ 

U 2 • R 2 ■ t/2' t/2 • T 3/2 and 


n+1/2 ET n+l/2 .„n+l/2 ...n+1/2 ..n+1/2 

P 3/2 * ET 3/2 ' EV 3/2 ' ** 3/2 ' ^ 3/2 


initially, and then sweeping the mesh as shown in Table V. 
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The quantities listed in (G.5) may be determined from equa- 
tions given previously, with the exception of the luminosity 
AL at the inner boundary. If this boundary is the center of 
the star, i.e., = 0 for all time n, then by symmetry 



G. 6 

for all n also. 

Another boundary condition must be imposed on diffusive 
energy transfer. Christy (1964) has proposed that this be 
accomplished by requiring that the surface boundary condition 
for the time-dependent problem be consistent with that for 
the time- independent diffusion equation. For a static, gray 
atmosphere, the solution of the equation of transfer for 
photons is 


T 4 = T fi 4 [t + q (T) ] 

where t is the optical depth and q(T) is a slowly- varying 
function. The diffusion approximation for the same problem 
gives 


T 4 - | T e 4 [T + c] 


where c is some constant . If c = 2/3 then 


»(T 4 > 




Surface 


3 4 

B — T 

4 e 


3 m 4 
2 T 


Surface 



If the effective temperature T q of the surface is known, 
the problem is determined. This may be expressed r.s 

- 4 !| 

' Surface 


5 (log T ) 
dT 


3 

2 


but to apply this expression it is necessary to know where 
the surface is. 

In order to avoid prejudging the calculation, a differ- 
ent approach was taken to determine the surface boundary 
condition. 

When the mean free paths per zone reached a certain small 
fraction X, a simple energy transfer calculation was made. 
The incident flux upon zone j+1/2 was obtained from the 
luminosity at boundary j while the opacity of zone j+1/2 
was determined by the temperature of this flux, i.e., the 
temperature of zone j-1/2. This gave the energy deposited 
while the neutrino loss rates discussed in chapter III 
gave the energy emitted by the zone j+1/2. Using (A. 13) 
with the change 

s n+l/2 _ AE n+l/2 
®j+l/2 j+1/2 

where ne ** ener 9y deposited by neutrinos, the 

temperature of zone j+1/2 in the emission surface was de- 
termined. In practice the transition from opaque to trans- 
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parent was so abrupt that the calculation was not sensitive 
to any sort of reasonable boundary condition of either of 
the types just mentioned. Outside the emission surface the 
uncoupled hydrodynamic scheme was used, without any diffusive 
energy transfer. 

H. Stability of the Difference Equations With Radiative 
Transfer . Richtmyer (1957) has discussed the stability of 
finite difference approximations to the diffusion equation 
in some detail. The discussion in this section is therefore 
limited to those aspects of stability of immediate interest. 

A complete treatment of the stability of a nonlinear diffu- 
sion equation coupled with the equations of hydrodynamics 
would be extremely complex. It appears that in practice 
the restrictions necessary for a linear, uncoupled diffusion 
problem can be suitably generalized for more complex systems. 

The simple form of the diffusion equation is 


be 

at 



in one dimension, where C is the concentration of whatever 
is diffusing, t the time, x the spatial coordinate and 0 
the diffusion coefficient. Perhaps one of the simplest 
difference approximations is 


? 
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c n+1 

_J 



At 


n+1/2 


D n+1/2 

3 



2C 11 + C 


/A 2 A n 

(Ax 


n 

2zl 


where the subscripts and superscripts have the same meaning 
as before, and 


<4x2) j = (ix j + l/2 + iX "-l/2 )2/4 


For stability it is necessary that 


H.l 


2DAt 

(Ax) 2 


for all j and n. This expression may be used to determine 
the time step At n+2/ ^ 2 at the next epoch. In particular, for 


the coupled problem, (H.l) gives 

n+1/2 n+1/2 ,n+l 


At 


n +3/2 . (R 1 


-J£l 


2 ( V " + /^) 2 ( T " : / 2 ) 3 


x const. 


'j-1/2 


j-1/2 ' 


H. 2 


which worked quite well when the minimum value for j = 1 , . . . 
JJ was taken. 


* See Richtmyer 1957, chapters 1 and 6 


7 


escape velocity for 
ejected matter 


"velocity of ejected 
matter" 


mass of remnant core 


radius of core 
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TABLE I 

Condensed model 
8 

~1.4 x 10 cm/sec 

9 

~7 x 10 cm /sec 

1.2 M 

o 

6 

1.6 x 10 cm 


kinetic energy of ^2 

ejected mass >4 x 10 ergs 


energy of emitted ^2 

electron neutrinos ~7 x 10 ergs 



? 
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TABLE II 


total light emitted 


intrinsic maximum 
visual magnitude 


mass ejected 


velocity of ejection 


« kinetic energy 


-17.5 


> 5 




g 

7 x 10 cm/sec 


>10 


52 


kinetic energy of 
■ ejected mass 


ergs 



S3 




TAB 

L E III 



(-* 


c~ & ^ 

r- h \ * 

3 l 

K 

decay./ 

V pair-annih.y 

V tablesy 

360 

5 

1rt 2 4 
x 10 

3 x 10 25 

1.4 x 10 26 

240 

3.7 

x 10 22 

4 x 10 23 

8 x 10 23 

120 

1.3 

x 10 21 

5 x 10 18 

2.6 x 10^ 8 


* These values are derived from the Chiu (1961) tables of 
electron pair- annihilation neutrino rates. 
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360 

240 

12G 


TABLE IV 


T 


10 


1.4 x 10 


10 


sec 


sec 


sec 



I 
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TABLE V 


quantity 

tiiae epoch 

space point 

U 

n+1/2 

j + 1 

R 

n+1 

j + 1 

V 

n+1 

j + 1/2 


n+1/2 

j +1/2 

T 

n+1/2 

j +1/2 

P 

n+1/2 

j+1/2 

ET 

n+1/2 

j + 1/2 

EV 

n+1/2 

j + 1/2 

AK • 

n+1/2 

j+1/2 

AL 

n+1/2 

j 

T 

n+1 

j-1/2 

P 

n+1 

j-1/2 
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Figure Captions 

1. Zero temperature equation at State. 

2. Comparison of the structure of isothermal-core 
and polytrope of index 3 models in the tempera- 
ture - density plane. 

3. Comparison of the structure of isothermal-core 
and polytrope of index 3 models in the radius - 
density plane. 

4. Evolutionary historyof M^_ = 1.5 zone of 10 M 0 
polytrope of index 3, for three different treat- 
ments of neutrino energy transfer. 

5. Evolutionary history of two representative mass 
zones of i 10 polytrope of index 3 initial 
model, with neutrino energy transfer treated in 
the diffusion approximation. 

6. Comparison with the calculations of Colgate and 
White. 

7. Temperature- density history of two representative 
zones of the isothermal-core initial model, with 


energy transfer by neutrinos in the diffusion 
approximation . 


Table Captions 


Table 

Table 

Table 


Table 

Table 


I. Comparison of centrally-condensed and poly tropic 
models. 

II. Observed characteristics of type II supernovae. 

III. Estimates of energy loss in ergs/gm/sec due to 

muon- type neutrinos formed by muon-decay and 

12 

muon pair-annihilation (density is p = 10 
gm/cm 3 ) . 

IV. Relaxation time for cooling due to muon- type 

neutrino escape at high temperature (density 
12 3 

is p = 10 gm/cm ) . 

V. Space- time points at which quantities appearing 
in the coupled difference equations of hydro- 
dynamics and diffusive energy transfer are 


evaluated 



